When can we detect that a P-f inite sequence is positive? 
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ABSTRACT 

We consider two algorithms which can be used for proving 
positivity of sequences that are defined by a hnear recurrence 
equation with polynomial coefficients (P-finite sequences). 
Both algorithms have in common that while they do succeed 
on a great many examples, there is no guarantee for them to 
terminate, and they do in fact not terminate for every input. 
For some restricted classes of P-finite recurrence equations 
of order up to three we provide a priori criteria that assert 
the termination of the algorithms. 

Categories and Subject Descriptors 

1.1.2 [Computing Methodologies]: Symbolic and Alge- 
braic Manipulation — Algorithms; G.2.1 [Discrete Mathe- 
matics]: Combinatorics — Recurrences and difference equa- 
tions 

General Terms 

Algorithms 

Keywords 

P-finite Sequences, Positivity, Cylindrical decomposition 

1. INTRODUCTION 

Inequalities for special functions are a serious challenge, 
both from the traditional paper-and-pencil point of view, 
but also (and in particular) for computer algebra. In con- 
trast to the vast number of algorithms for dealing with iden- 
tities, almost no algorithms are available for inequalities. 
Already for the very restricted class of sequences satisfy- 
ing linear recurrence equations with constant coefficients 
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(C-finite sequences), the positivity problem leads to hard 
number theoretic questions to which no solutions are known 
today, see [8, 10] and the references given there for the cur- 
rent state of the struggle. 

Still, inequalities are not entirely hopeless. For example, 
Mezzarobba and Salvy have recently given an algorithm for 
effectively computing tight upper bounds for sequences de- 
fined by linear recurrence equations with polynomial coeffi- 
cients (P-finite sequences) [16]. Five years ago, Gerhold and 
Kauers [11] proposed a method applicable to inequalities 
concerning quantities that satisfy recurrence equations of a 
very general type. Their method consists of constructing a 
sequence of polynomial sufficient conditions that would im- 
ply the non-polynomial inequality under consideration. If 
one of the conditions in the sequence happens to be true 
(which can be detected, e.g., with Cylindrical Algebraic De- 
composition [5, 6, 4, 2]), the method succeeds, otherwise 
it keeps on running forever. Simultaneously, the method 
searches for counterexamples and it will find one and termi- 
nate for every false inequality. 

Despite its simplicity, the method has proven quite suc- 
cessful in applications. Not only did it provide the first 
computer proofs of some special function inequalities from 
the literature [11, 12, 13, 14], but it even helped to resolve 
some open conjectures [1, 15, 14, 17]. At the same time, the 
method remains somewhat unsatisfactory from a computa- 
tional point of view, as it is not clear on which inequalities 
it succeeds and on which it doesn't. It would be interesting 
to have, at least for some restricted classes, some a priori 
criteria telling us whether the method (or some variation of 
it) will succeed or not. 

Our goal in this paper is to provide such criteria for two 
particular proving procedures (Algorithms 1 and 2 described 
below). We are far from being able to give a full answer to 
the question posed in the title, but we can identify some 
nontrivial portions of P-finite recurrence equations of fixed 
order on which termination of Algorithms 1 or 2 is guaran- 
teed. For first order equations, deciding positivity is trivial. 
For second order equations, we provide a result (Theorems 2 
and 3) that answers the question under a genericity assump- 
tion. For third order equations, we are able to identify the 
terminating cases of Algorithm 2 but only have partial re- 
sults for Algorithm 1 supplemented by empirical evidence 
supporting a conjecture concerning its terminating cases. 
An interesting aspect of our analysis is that algorithms for 
real quantifier elimination are not only used as a subroutine 
of Algorithms 1 and 2, but they are also contributing in an 
essential way to the proofs of our termination theorems. It 



is therefore possible — in principle — to extend our results to 
equations of order greater than three. Only the increasing 
time and memory requirements of the computations have 
prevented us from doing so. 

2. PRELIMINARIES 

A sequence /:N— >■-?(" :=]Rn(Q is called P- finite (or 
holonomic) if there exist polynomials po,. . . ,pr £ K[x\, not 
all zero, such that 



po{n)f{n) +pi(n)/(n + l) 



+ Pr{n)f{n + r) = 



for all n G N. Such an equation is called a (P-finite) recur- 
rence, and r is called its order. If pr("-) 7^ for all n € N, 
then the infinite sequence / is uniquely determined by the 
recurrence and r initial values /(O), /(I), • • • , f{r — 1). The 
assumption pr (n) 7^ for all n £ N can be adopted without 
loss of generality, because we can substitute g{n) — f{n + u) 
for some u larger than the greatest integer root of pr and 
then consider g instead of / and check nonnegativity of the 
finitely many terms /(O), /(I), • • • , f{u — 1) by inspection. 
We will do so. 



From now on, all recurrences are assumed to 
have a leading coefficient pr with no positive in- 
teger roots. 



A P-finite recurrence is called balanced if degpo ~ deg pr 
and deg Pi < degpo (* ~ 1, ■ • ■ , '')• The characteristic poly- 
nomial of a balanced recurrence is defined as 

^Cy{po{y) + pi{y)x -\- p2{y)x'^ -\ \-Pr{y)x'') G K[x]. 

Its roots Ai,...,Ar G C are called the eigenvalues of the 
recurrence. (The \i are not necessarily distinct.) Note that 
for a balanced recurrence, the characteristic polynomial has 
always degree r and it has never as a root. 

An eigenvalue A; is called dominant if |Aj| < |Ai| for all 
j = 1, . . . ,r. Dominant eigenvalues govern the asymptotics 
of the sequences defined by the recurrence [20, 9]. If there 
is a unique dominant eigenvalue Ai, then for we will usually 
have 

f{n) ~ c(n)A" (n -^ 00) 

where c is of subexponential growth in the sense that 

c{n-\-l) n- 



c{n) 



^ 1. 



There may be choices of initial values for which c{n) = 
for all n so that the asymptotics of / is not affected by Ai 
but by the next smaller eigenvalue(s). Whether this is the 
case or not can be hard to verify formally, but is usually 
easy to verify empirically. Some of our termination results 
apply only to this generic situation where initial values are 
chosen such as to actually exhibit the asymptotic behavior 
predicted by the dominant eigenvalue. 

Finally, if the dominant eigenvalue Ai is not real and pos- 
itive, then it is clear that / will be ultimately oscillating, 
and so fin) > cannot possibly be true for all n. This case 
can be sorted out trivially beforehand, and we may there- 
fore assume that the unique dominant eigenvalue is real and 
positive. In this case, the substitution g{n) = f{n)/X" turns 
the recurrence 

Po{n)f{n) + pi{n)f{n + 1) + ■ ■ ■ + Pr{n)f{n + r) = 



into 

po{n)g{n) -|-pi(n)A,,g(n + 1) + ■ 



+ Pr{n)Xlg{n-\-r) = 



whose dominant eigenvalue is 1. As g{n) > <=> .f{n) > 
0, it suffices to consider this case. 

3. INDUCTION BASED PROVING 
PROCEDURES 

3.1 The Original Version 

The approach of [11] is as follows. Suppose that /: N — ^ 
K is defined by a recurrence 

Po{n)f{n) +pi(n)f{n +!) + ■■ ■+pr(n)f{n + r) ^ 

and initial values /(O) = /o, /(I) = /i, . . . , /(r - 1) = /,.-i. 
We seek to prove f{n) > for all n £ N by induction: 

/(n) > A ■ • ■ A /(ti + r - 1) > ^ /(n + r) > 0. 

Because of the recurrence, this is equivalent to 

f{n) > A • ■ • A /(n + r - 1) > 



Po{n) 



fin) 



Pr-i{n) 



/(n + r-1) > 



Pr{n) Pr(n) 

For this to be true for all n G N, it is sufficient that the 
induction step formula 

y yo,yi,...,yr-i G IR V 2; G IR : 
(x > A yo > A • ■ • A yr-i > O) 



Pojx) 

Pr{x) 



ya 



iW 



Pr{x) 



J/r-l > 



is true, and this can be decided by a quantifier elimination 
algorithm. If it is true, the induction step is established and 
/ is nonnegative everywhere if and only if it is nonnegative 
for n — 0, . . . ,r — 1, which can be checked. 

In the unlucky case when the induction step formula is 
false, there is no immediate conclusion about / that could 
be drawn. In this case, refined induction step formulas 

f(n) > A ■ • ■ A /(ti -h £. - 1) > =^ /(ti -f £.) > 

for g > r are constructed. Using the recurrence, each term 
f{n -\- i) can be rewritten as a linear combination of f{n), 
. . . , f{n-\-r — l) with rational function coefficients, and using 
this rewriting, the refined induction step formula takes the 
form 

^{g) :=Vyo,j/i,...,y,-i G K, V x G IR : 

(a- > A yo > A ■ ■ • A yr-i > 

A qrfi{x)yo H h qr,r-i{x)yr-i > 

A qr+ifiix)yQ -\ + qr+i,r-i{x)yr-i > 



A qg-i.o{x)yo H + qg-i,r-i{x)yr^i > O) 

=> qgfi[x)yo H + qe,r~i{x)yr-i > 0, 

where the qi^j are some rational functions. 
The full method then reads as follows. 



Algorithm 1. 
Input: A P-finite recurrence of order r and a vector of initial 



values defining a sequence / : N —> A". 

Output: True if f{n) > for all n £ N, False if /(n) < for 

some n G N, possibly no output at all. 

1 for n = to r — 1 do 

2 if f{n) < then return False 

3 for n — r,r + l,r + 2,r + 3, . . . do 

4 if <l>(7i) then return True 

5 if f{n) < then return False 



return False 
else if 3 /^ > 



Example 1. Let /: N — s 

{2n + 13) fin + 3) - 
+ (Sn + 20) f{n H 
/(0) = /(l)=/(2) 



K be defined by 

(5n + 22) f{n + 2) 
1) - (2n + 7) fin) 



0, 



We use Algorithm 1 to show that fin) > for all n £ N. 

Since /(O), /(I), /(2) > 0, we enter the loop in line 3. For 
n = 3, we have 

<l>(n) =V j/o,yi,y2 V x e R: 

(a: > A yo > A yi > A y2 > 0) 



mhy^ 



ta+20 , 5^+22 > 

21 + 13 yl ^ 2a; + 13y2 ^ "• 



T/iis is /aZse, 6tii since /(3) = 9/13 > (checked in line 5), 
we continue. 

The formula $(4) is too lengthy to be reproduced here ex- 
plicitly, and it is also false. Yet /(4) = 61/195 > 0, so we 
proceed to consider the even lengthier formula $(5), which 
turns out to be true. At this point the algorithm terminates 
with output True. 

3.2 A Variation 

In cases where Algorithm 1 does not terminate, it is some- 
times possible to prove inductively the stronger statement 
that fin) is increasing, viz. that /(n + 1) > fin) for all 
n > 0. While this is obviously a sufficient condition for 
fin) > for all n, there are of course sequences / which 
are non-negative but not increasing. For such cases, a good 
strategy is to prove that /i~"/(n) is increasing, for some 
suitably chosen constant fi > 0. The choice of ^ is criti- 
cal in two respects: it must be small enough to assure that 
/x~"/(n) actually is increasing, and it must be big enough to 
allow for an inductive proof. The following algorithm proves 
positivity of a P-finite sequence / by searching for a fi that 
meets both criteria. 



Algorithm 2. 
Input: A P-finite recurrence of order r and a vector of initial 
values defining a sequence / : N — >■ A". 

Output: True if /(n) > for aU n £ N, False if /(n) < for 
some n G N, possibly no output at all. 

1 Determine a quantifier free formula ^i^, n) equivalent 
to 

V yo,...,yr-i V a: > C : 

f j/o > A yi > nyo A ■ • ■ A i/r-i > A*yr-2 j 



Po(x) 



tVo 



_ Pr-lix) 
Prix)"" Prix) 

2 for n = 0,1,2,3,... do 

3 if fin) < then 



yr~l > f^Vr-l 



A- 



$(n,/i)A/(n-H)>M/(n) 
A/(n + r-l) > /i/(n + r - 2) then 



6 return True 

Theorem 1. Algorithm 2 is correct. 

Proof. Correctness is obvious whenever the algorithm 
returns False, because this happens only when an explicit 
point n with fin) < has been found. Suppose now that the 
algorithm returns True at the nth iteration of the for loop. 
Then /(fc) > for k = 0, . . . , n, otherwise the algorithm 
would have terminated in an earlier iteration with output 
False. The condition in line 5 inductively implies 

3/i>0Vfc>n: fik + 1) > fifik). 

Since /i > and fin) > 0, this inductively implies /(fc) > 
also for all fc > n. D 

Example 2. Let /: N — > A' be defined by 

in + 3) fin + 3) - (5n + 13) fin + 2) 

+ (5n + 12) fin + 1) - in + 2) fin) = 0, 
/(O) = 1, /(I) = 1/4, /(2) = 1/10. 

Algorithm 1 does not seem to terminate for this sequence. 
But Algorithm 2 succeeds. 

Step 1 produces the quantifier free formula 



,.„^5-^/5^ ^ ysg^ + 22e + 25 -f 5g + 13 
C>OA^^<M< ^^^^^^ , 

which we denote $(^,/i). In the iteration of the for loop, we 
get: 

For n = 0, since /(O) = 1 > 0, uie check whether 

0>0A^^<M<3Ai>MA^>J 

is satisfiable. As it is not, we proceed. 

Also n — 1 and n — 2, we have fin) > but there is no 
jj. > with 

$(n, ^) A fin + 1) > ^/(n) A fin + 2) > fifin + 1). 

Then for n — 3, since /(3) > 0, we check whether 

3>0A^^^<M<^^ + ^^AlI>iiAgiI>lZ^ 
2 - ^ - 12 80 - 10 400 - 80 

is satisfiable. As it is satisfiable (e.g., by fi — 2) the algo- 
rithm terminates with output True. 

The two strategies employed in Algorithms 1 and 2 in the 
case where a direct proof of the induction step formula fails 
(prolonging the induction hypothesis in Algorithm 1 ver- 
sus multiplying with a positivity preserving exponential in 
Algorithm 2) are independent of each other. It is possible 
to merge both strategies into a single strategy that simul- 
taneously prolongs the induction hypothesis and inserts a 
positivity preserving exponential. An algorithm based on 
such a combined strategy is easily written down, but turns 
out to be computationally quite expensive on examples. It 
would be interesting to carry out the termination analysis 
given below for the combined algorithm, but the quantifier 
elimination problems arising in this analysis seem currently 
too hard to be carried out for the combined algorithm. 



4. TERMINATING CASES 

Both algorithms given in the previous section may fail to 
terminate. Our goal now is to identify classes of P-finite re- 
currence equations for which termination can be guaranteed 
a priori. 

4.1 Order One 

This case is rather simple and included here merely for 
the sake of completeness. If /: N — )> iv satisfies 

po{n) f{n)+pi{n)f {71 +1)=0, 

then /(n) > for all n G N if and only if /(O) > and 
—po(n)/p\(n) > for all n G N. Since sign changes of 
—po{n)/pi{n) can occur only at the real roots of po or pi, 
the only thing we need to do is to find an upper bound 
no G IR for the real roots (this can be done) , and then check 
whether —po{n)/pi{n) > for n = 0, 1, 2, . . . ,no + 1. 

Example 3. Consider f : ¥i -^ K defined via 

(3n - 16)/(n) - (3n - 17)/(n + 1) = 0, /(O) = 1. 

The roots ofpo,pi are 16/3 and 17/3, respectively, and they 
are both less than no — 6, for instance. Therefore, since 
fin) > for ri — 0, ... ,6, we can conclude that f{n) > 

for all n G N. 

4.2 Order Two 

We now turn to sequences / : N — >■ _R' which are defined 
by a balanced P-finite recurrence of second order, 

P2(n)/(n H- 2) - pi(n)/(n + 1) - po{n)f{n) = 0. 

We assume (without loss of generality) that 1 is a dominant 
eigenvalue of this recurrence and let u £ K with |m| < 1 be 
such that 

{x ~ \)(x — u) = X — (u + V)x — {~u) 

is the characteristic polynomial of the recurrence. The ques- 
tion is whether Algorithm 1 and Algorithm 2 succeed in 
proving that fin) > for all n. (If this is actually the case; 
if it is not, then both algorithms will obviously succeed in 
finding a counterexample.) We will show that termination of 
Algorithm 1 depends on the sign of u whereas Algorithm 2 
(generically) terminates for all u. 

Theorem 2. Ifu G (—1,0), then Algorithm 1 terminates. 
Proof. Rewrite the recurrence in the form 

/(n + 2) = ^/(n + l) + ^/(n). 
P2{n) P2{n) 

Since the characteristic polynomial is 

(a; — l){x — u) = X — {u+ l)x — {^u), 

we have 

Pl(n) n^oo , 1 .^ n ] Po{n) n^oo „ 

; : > u + 1 > and ; ; > -u > 0. 

P2(n) P2(n) 

Therefore, 

3n„GNVn>no:^>0A^>0, 
P2{n) P2{n) 

where we may safely regard n as ranging not only over the 
integers but over all reals except the roots of P2- We show 



that Algorithm 1 terminates after at most no iterations. The 
previous formula implies 

V yo , yi V n > no : (yo > A j/i > O) 

P2(n) P2(n) 

Substituting n i— >■ n — no leads to 

Vyo,yi Vn>0: (i/o> OAyi >0) 

^ Po(n + no) pi(n-l-no) 

:y" "I — 7 — ; \ 2/1 ^ 0- 



P2(n + no)' 



P2{n + no) 



As the variables j/o, yi range over all reals, the latter formula 
will remain true if we apply a substitution 

yo 1-^ ro(n)yo + ri(n)j/i, 
2/1 !-)> r2{n)yo + r3{n)yi, 

where ro, . . . , rs are some rational functions in n. This gives 
V yo,yi V n > : {ro{n)yo + ri{n)yi > 

A r-2(n)j/o + r-i{n)yx > O) 
Po(n-|-no)ro(n) -|-pi(n -|- no)r2(n) 



+ 



/ , ^ "2/0 

P2(n-|-no) 

po(n + no)ri(n) +pi(n-f no)r3(n) 

P2(n + no) 



yi>o. 



We are free to further modify this formula, without harming 
its truth, by imposing arbitrary additional conditions on the 
left hand side of the implication. 
By choosing ro, ri, r2, rs such that 

/(n + no) = ro(n)/(n) + ri(n)/(n + 1) 
/(n + no + 1) = r2(n)/(n) + r3(n)/(n + 1) 

and adding constraints qi,o{n)yo H- (?i,i(n)j/i > encoding 
fin -\- i) > (i = 0, . . . , no — 1) on the hypothesis part, we 
obtain precisely the formula $(no) as defined in Section 3.1 
and used in Algorithm 1. As the formula is true, the al- 
gorithm terminates in the noth iteration (or earlier), as we 
wanted to show. D 

Remark 1. Algorithm 1 fails to terminate for positive u. 
To see this, consider the C-finite recurrence 

fin + 2)-(u+ l)/(n + 1) + ufin) = 

for some u G (0, 1). If Algorithm 1 applied to this recurrence 
terminated in the noth iteration, for some no > 0, then the 
truth of $(no) implies that no solution f:¥S^>-Kof the 
recurrence can have no consecutive nonnegative terms fol- 
lowed by a negative term. (So that, if no consecutive terms 
are found nonnegative, all subsequent terms must be non- 
negative as well.) 

To see that no such no can exist for the C-finite recur- 
rence above, it is sufficient to construct for every no > 
a solution which contains a run of exactly no nonnegative 
terms. The general solution of the recurrence is co -\- cim" 
for some constants co,ci £ K . It is easily checked that the 
choice Co = —1, ci = u"""^ has the desired property. 

The argument extends, at least for generic initial values, 
to P-finite balanced recurrence equations, using the fact that 
the recurrence admits two solutions /i , /2 : N — s> K with 
/i(^+l)//iW ^^^^^ 1 and/2(n + l)//2(n) -^^^^ u. 



Theorem 3. If u e (—1, 1) \ {0}, then Algorithm 2 ter- 
minates for generic initial values. 

Proof. Consider the set Da C ]R^ consisting of all points 
(co , ci , /i) satisfying 

0<At<lA^<ci<2A ^(^ - ci) < CO < 1 

and the set 

D2 ■■= { (co, ci) e R^ : < ci < 2 A -\cl < co < 1 }. 

It can be shown by CAD computations that 

V (co,ci) G D2 3 mG (0,1) : (co,ci,^) G Dz (1) 

and that 

V (co,ci,^) G -D3 V yo,yi G K, : 

(yo > A j/i > ^yo) => coyo + ciyi > ^j/i. (2) 

Since i(u + l)^ > m for all u G ( — 1, 1), the set D2 contains in 
particular the point (— u, k+1) where u is from the statement 
of the theorem. Because of (1), there exists /i G (0, 1) with 
{—u,u + l,fi) G D3. Since D3 is open, there exists e > 
such that 

U ;= (— u ~ e,—u + e) x (u + 1 — e, u + 1 + e) x {^} C D3. 

Because of 



Po{n) n^oo , Pi(n) n- 

> —u and 



P2(n) 



->■ w + 1, 



,/i eUCDs 



P2{n) 

there exists ^ G N such that 

/ PoW Pi{n) 
\ P2(n)' p2(n)' 

for all n > ^. Together with (2), this implies 

3/iG(0, l)3CGNVn>CV yo, yi G IR : 

/ \ Voin) Po{n) 

[yo>OAyi> p-yo) => ,_^ yo + ,_^ yi > pyi- 



P2(n) 



P2(n) 



Therefore, the set 

C:={(e,M)e(0,(»)x(0,l):$(e,^)} 

with $(^, /i) as used in Algorithm 2 is not empty. 

Fix some point (^, fj,) G C. Then it is immediate by the 
defining formula that also (^', p) £ C for every ^' > 5, so 

(C,«))x{M}ca 

Let now /: N — > iiT be the sequence to which Algorithm 2 
is applied. Then, because the eigenvalues of its defining 
recurrence are 1 and u and we have |m| < 1 and we assume 
generic initial values, we have 

/(n+1) 



fin) 



^ 1. 



Since fi < 1, this implies the existence of an index m G N 
such that f{n) > for all n > m and f{n + l)//(n) > /x for 
all n > m, so that we get f{n + 1) > fJ-f{n) for all n > m. 

It follows that the algorithm terminates no later than at 
iteration max(m,^). D 

Remark 2. The defining inequalities of D3 m the pre- 
ceding proof were found by quantifier elimination applied to 
formula (2) with the first quantifier dropped. This computa- 
tion as well as the CAD computations referred to m the proof 



were performed with Mathematica's built-in implementation 
of CAD [18, 19]. The computation time is negligible for all 
of them and we are sure that other implementations [7; 3, 
etc.] would have no problem with them either. 

Example 4. The restriction to generic initial values in 
Theorem 3 is essential: let f : M -^ K be defined via 

{n + 3ff{n + 2) - i(n + 2)(3n + ll)/(n + 1) 
+ i(n + 4)(n + l)/(n) =0 

and /(O) = 1, /(I) = 1/4. Then we have f{n) = 2-"/(n + l) 
for all n £M and so in particular f{n) > (n G Nj. 
Algorithm 2 finds 



$(C,/i) = i<M<lAS> 



17- 12^ - V25 - 16m 
2/i-3 



in Step 1 and continues by searching for an index n with 
f{n + 1) > ^f{n). As no such index exist, the search con- 
tinues forever. 

The general solution of the defining recurrence for f is 



Co + Cl 



2"" 
n+ 1 



and we will have co 7^ for a generic choice of initial values. 
In these cases, the solution converges to Co and therefore 
eventually reaches an index no with a term that is greater 
than half its predecessor. 

4.3 Order Three 

Consider now sequences / : N — >■ i^ defined by a balanced 
P-finite recurrence of third order, 

P3(n)/(n+3)-p2(n)/(n+2)-pi(n)/(n+l)-po(n)/(n) = 0. 

Again, we assume without loss of generality that 1 is a dom- 
inant eigenvalue and we let u,v £ K he such that 

(x — l)(a:: -{- ux -'r v) — x — (1 — u)x — {u — v)x — v 

is the characteristic polynomial of the recurrence under con- 
sideration. The condition that the two roots of the quadratic 
factor belong to the interior of the complex unit disc trans- 
lates into the condition 

\u\~Kv <1 

for the coefficients of the polynomial. The points {u, v) G K^ 
satisfying this condition form the interior of the triangle with 
corners at (-2, 1), (2, 1), (0, -1): 




-1012 
Just for the sake of orientation: the polynomial x'^ +ux-\-v 
has two complex conjugate roots in region A, two positive 
real roots in region B, two negative real roots in region C, 
and a positive as well as a negative root in region D. 



We want to identify regions of the triangle corresponding 
to recurrence equations on which Algorithms 1 and 2 termi- 
nate. Only for Algorithm 2 we have a satisfactory result, so 
let us consider this case first. 

Theorem 4. //|u| — 1 < v < 1 andAv < (u+1)^ andu < 
1, then Algorithm 2 terminates for generic initial values. 

Proof. Consider the set D4 C ]R* consisting of all points 
(co , ci , C2 , ^) satisfying 

0</i<lA^<C2A /^(/i — C2) < ci A /i^ — C2/i^ — cifi < Co 

and the set D3 C ]R consisting of all points (co,ci,C2) with 

(0 < C2 < 2 A -ic2 < ci < min(3 - 2c2,cl) 

A 2ci + 9ciC2 + 27co + 2{cl + Saf^'^ > O) 

V (0 < C2 < 1 A ci > C2 A Co + C1C2 > 0) . 

The following facts can be verified by CAD: 

• V (co,ci,C2) G D3 3 ^ G (0, 1) : (co,ci,C2,m) £ D4 

• V (co,ci,c2,/i) G -D4 V yo,yi,y2 G IR : 

(yo > A yi > ^yo A 1/2 > ^iJ/i) 

^=> Coj/o + ClJ/l + C2y2 > A*J/2 

• V u, -u G R : (|u| - 1 < u < 1 A 4u < (u + 1)^ A u < 1) 

^=> {v,u — u, 1 — u) G -Dg- 

Consequently, for ti, v from the statement of the theorem 
there exists n G (0, 1) such that {v,u — v,l — u,fi) G D4. 
Since D4 is open, there exists e > such that 

U := {v — e,v + e) x {u — v — e,u — v + e) 

X {I — u — £.1 — u + e) X {^} C D4. 



Using 

PoM n->oo Pl(w) n->cx> 



P3(n) 



P3(ri) 



P2{n) r^ 

pain) 



-^ I — u, 



the rest of the proof is fully analogous to the proof of The- 
orem 3. n 

The set of points (u, v) for which Theorem 4 asserts ter- 
mination of Algorithm 2 is the shaded area in the figure 
below. 




-2 -1 

The truth of the formula 

V-uG (-1,1) V-uG (1^1-1,1) : 
[3 /i > V yo,yi,y2 ■ 

(yo > A yi > nyo A y2 > ^J.yl) 
^=> vyo + (u — v)yi + (1 — u)y2 > fJ,y2] 
=> u<l A4v < {u + if 



(as confirmed, once again, by a CAD computation) asserts 
that Theorem 4 is sharp. 

We are not able to provide a sharp result for the ter- 
minating region of Algorithm 1. If we proceed to reason 
as in the proof of Theorem 2, we obtain termination for 
{u,v) restricted to the (open) triangle with vertices (0,0), 
(1,0), (1,1), essentially because of 

Vug (0, 1) V«G (0,m) Vyo,yi,y2 : 

(yo > A yi > A y2 > 0) 

=^ vyo + {u — i')yi + (1 — u)y2 > 

and the convergences 

Poin) n^oo Plin) n^oo P2{n) 



p-i{n) 



P'iin) 



ps{n) 



^ 1 



But this is not the entire terminating region. A larger por- 
tion of the terminating region can be identified by starting 
out with a formula corresponding to an induction hypothesis 
of length four. As the formula 

V yo, yi, y2 : (t/o > A yi > A y2 > 

A vyo + (m — v)y-i + (1 — u)y2 > O) 

^=> (1 — u)vyo + u(l — u + v)yi + (1 — it -f u — v)y2 > 

is true for all {u,v) with 

u<lAv>QM-u + u^-v>G 

A {u > QV u^ - V - uv + v^ < Q) , 

and as we have 

Po{n)p2{n-\- 1) 



^> (1 — u)v, 

-->■ u(l — It + «), 



-> 4 — u + u 



P3{n)p3{n + 1) 
Pi {n)p2 (w -I- 1) -I- po (n + l)p3 (n) ^ 

P3{n)p3{n + 1) 
P2 {n)p2 (n -I- 1) -I- pi (n -|- l)p3 (n) ^ 

P3{n)p3{n + 1) 

Algorithm 1 also terminates for all (it, v) satisfying the con- 
ditions stated above. 

Starting out with a formula corresponding to a hypothesis 
of length five leads to a portion of the termination region 
whose description can be computed in a reasonable amount 
of time, but which is already too big to be reproduced here. 
For longer induction hypotheses, the computational effort 
for doing quantifier elimination becomes prohibitive. But 
it is still possible to determine experimentally the regions 
obtained by taking a particular length q of the induction 
hypothesis taken as the starting point of the termination 
proof. The empiric results for induction hypotheses of length 
up to 10 are as follows (the numbers indicate the length of 
the induction hypothesis): 




The picture suggests the following characterization for the 
full region of termination. 

Conjecture 1. //|uj — 1 < u < 1 and 

{u> lAv > 0)V 4v > {u+ if, 
then Algorithm 1 terminates. 

The conjecture is equivalent to saying that Algorithm 1 
terminates if i^ + ita; + u has no positive root. If the conjec- 
ture is true, then about 96.35% of the area of the triangle 
are covered by one of the two algorithms we considered. 
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